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INVESTIGATION OF THE LAMINAR AERODYNAMIC HEAT-TRANSFER CHARACTERISTICS 
OF A HEMISPHERE-CYLINDER IN THE LANGLEY 11-INCH HYPERSONIC TUNNEL 

AT A MACH NUMBER OF 6.8 f 

By Davis H. Crawford and William D. McCauley 


SUMMARY 

A program to investigate the aerodynamic heat transfer of a 
nonisothermal hemisphere-cylinder has been conducted in the 
Langley 11 -inch hypersonic tunnel at a Mach number of 6.8 
and a Reynolds number from approximately 0.14- X 10 6 to 
1.06 xlO 6 based on diameter and free-str earn conditions. The 
experimental heat-transfer coefficients were slightly less over the 
whole body than those predicted by the theory of Stine and 
Wanlass ( NACA Technical Note 3344) for an isothermal 
surface. For stations within 4 6° of the stagnation point the 
heat-transfer coefficients could be correlated by a single relation 
between local Stanton number and local Reynolds number. 

Pitot pressure profiles taken at a Mach number of 6.8 on a 
hemisphere-cylinder have verified that the local Mach number or 
velocity outside the boundary layer required in the theories may 
be computed from the surface pressures by using isentropic flow 
relations and conditions immediately behind a normal shock. 
The experimental pressure distribution at a Mach number of 
6.8 is closely predicted by the modified Newtonian theory. The 
velocity gradients calculated at the stagnation point by using 
the modified Newtonian theory vary with Mach number and 
are in good agreement with those obtained from measured pres- 
sures for Mach numbers from 1.2 to 6.8. 

At the stagnation point the theory of Sibulkin, in which the 
diameter and conditions behind the normal shock were used, 
was in good agreement with the experiment when the velocity 
gradient at the stagnation point appropriate to the free-stream 
Mach number was used. 

INTRODUCTION 

One of the earliest experimental investigations in connec- 
tion with heat transfer on a hemisphere-cylinder at super- 
sonic speeds was carried out by Eber (ref. 1) who obtained 
the “brake” temperature distribution along the surface. 
These brake or effective recovery temperatures decreased 
with increasing distance from the stagnation point on the 
hemisphere and were influenced by the free-stream Mach 
number. Ivorobkin (ref. 2) and Stine and Wanlass (ref. 3) 
show the same type of distributions and, in addition, demon- 
strate that for the hemisphere-cylinder with a laminar 
boundary layer the recovery factor r) r may be expressed by 
the square root of the Prandtl number based on the local 
conditions just outside the boundary layer. 


Stabler and Nielsen (ref. 4) compared the average heat 
transfer (Nusselt number) with Reynolds number for differ- 
ent- Mach numbers and found that, in their experiments, 
data at subsonic and low supersonic Mach numbers could be 
correlated when the supersonic data were based on conditions 
behind the normal shock. Experimental investigations of 
the local heat transfer on a hemisphere-cylinder have been 
made by Korobkin (refs. 2 and 5) and Stine and Wanlass 
(ref. 3). Stine and Wanlass noticed that, as the free-stream 
Mach number was increased, the pressure distributions were 
tending toward a common curve and also suggested that the 
heat-transfer results they obtained at a Mach number of 
1 .97 could be representative of those obtained on similar 
isothermal surfaces with laminar boundary layer for Mach 
numbers greater than 1.97 and of temperatures less than 
that of dissociation. Although this tends to be true, a Mach 
number of 1.97 is probably too low. 

Some theoretical methods for obtaining local heat transfer 
on a hemisphere-cylinder at supersonic speeds have as their 
bases the solutions of incompressible flow on isothermal sur- 
faces. An exact incompressible solution utilizing a method 
applied by Squire (presented in ref. 6) to a cylinder was given 
by Sibulkin in reference 7 for the heat transfer at the stagna- 
tion point. Sibulkin suggested that this solution may be 
applicable at supersonic speeds when conditions behind the 
normal shock are used. An approximate solution was 
derived for incompressible flow over a hemisphere with an 
isothermal surface by Sibulkin (presented in ref. 2) by 
assuming that the velocity and temperature profiles were 
similar and of parabolic shape. Korobkin (ref. 2) compared 
these theoretical results with his data at a Mach number of 
2.80 and found that, when the incompressible-flow solution 
was adjusted so as to pass through Sibulkin’s stagnation- 
point value, it roughly approximated the data. However, 
this result does not agree with the theory by Stine and Wan- 
lass (ref. 3) who used a different approach. By using the 
Mangier transformation, the Stewartson transformation, 
and thermal solutions to the Falkner-Skan wedge-flow 
problem, they evaluated the heat-transfer rate in axisvm- 
metric flow for an isothermal surface in terms of the known 
heat-transfer rate in an approximately equivalent two- 
dimensional flow. 

The purpose of this investigation was to extend the range 


tSuperscdes NACA Technical Note 370G by Davis H. Crawford and William D. McCauley, 1966. 
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of the previous investigations for heat transfer on hemisphere- 
cylinders to a Mach number of 6.8. This investigation 
utilized data from transient noniso thermal temperature dis- 
tributions since, at the high temperature levels involved, an 
isothermal surface was difficult to obtain. 
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SYMBOLS 


speed of sound 
semiliemispherical arc length 
specific heat of model wall material 
specific heat of air at constant pressure 
specific heat of air at constant volume 
diameter 


local aerodynamic heat-transfer coefficient, T ~ 

-L r J- s 

thermal conductivity for air 
thermal conductivity for model 
a characteristic length 
Mach number, u/a 
Nusselt number, hL\k 
Prandtl number, c P nlk 
pressure 


heat flow to model per unit 
Reynolds number, 


area, — k 



radius 

Stanton number based on x and local free- 
stream conditions just outside boundary layer, 
A r s t = N Nu /RN Pt = h/pucp 
absolute temperature 
velocity 

distance along body generatrix from stagnation 
point 

space coordinate normal to body 


velocity gradient along surface, 


bu 

bx 


emissivity 

. * T r -T l 
recovery factor, — =• 

angle from stagnation point on sphere 

distance along cylinder from juncture 

absolute viscosity 

kinematic viscosity 

angle about axis of model 

mass density of air 

mass density of model shell material 

Stephen-Boltzmann constant 

time 


Subscripts: 

a condition of model's surroundings 

D parameter based on model diameter as a char- 

acteristic length 

h pitot pressure behind bow shock 

i not. corrected for heat conduction along model skin 

in local condition on inside surface of model shell 

L parameter using L as characteristic length 

r recovery, surface temperature 

s local conditions on outside surface of model 


t stagnation conditions 

x parameter based on static conditions just outside 

boundary layer and with x as characteristic 
length 

<t free-steam conditions immediately behind normal 

shock 

l local free-stream conditions just outside boundary 

layer 

oo undisturbed free-stream conditions 

APPARATUS 

TUNNEL 

The present investigation was conducted in the Langley 
11-inch hypersonic tunnel described in references 8 and 9. 
Air is stored at 50 atmospheres pressure and is released 
through an adjustable pressure regulating valve and a new 
direct air heater with tube resistance elements of nickel- 
chromium alloy (replacing the storage heater described in 
refs. 8 and 9) to the settling chamber and nozzle. A two- 
dimensional nozzle constructed of Invar and designed for a 
Mach number of 7 was used for this investigation. Invar 
was used in the construction of this nozzle in order to alleviate 
the deflection of the first minimum, which occurred in the 
steel nozzle because of differential heating of the nozzle 
blocks. Preliminary calibrations have indicated a Mach 
number of 6.86 ± 0.04 at 30 atmospheres stagnation pressure 
in the central core of uniform flow which measures about 
6 y 2 inches in the vertical direction by about 6 inches in the 
horizontal direction. 


MODELS 

Four models were used in this investigation. A cross- 
sectional view showing construction details typical of all 
models is presented in figure 1. Photographs of the separate 
models are shown in figure 2. Two models, one of 1.170-inch 
diameter and the other of 3.025-inch diameter, were instru- 
mented to provide transient temperature data, a third model 
was instrumented to provide pressure data, and a fourth 
model, temperature-recovery data. Orifices and thermo- 
couples were installed on the individual models at stations 
along the body generatrix as indicated in table I. The heat- 
transfer and pressure-distribution models were constructed 
of SAE 1020 carbon steel and the temperature-recovery model 
was constructed by electroplating nickel on a mandrel. The 
dimensions of the heat-transfer models were carefully meas- 
ured before the parts were assembled with silver solder. 

The small heat-transfer model was 1.170 inches in diameter, 
5.32 inches in length, and varied from 0.030 inch to 0.035 inch 
in wall thickness. Temperature measurements on this model 
were taken from eight chromel alumel thermocouples of 
No. 30 gage wire, five of these being located at stations from 
the stagnation point to the hemisphere-cylinder juncture. 
(See table I.) 

The large heat-transfer model was 3.025 inches in diameter, 
10.5 inches in length, and varied from 0.098 inch to 0.101 inch 
in wall thickness. It contained 25 thermocouples, silver sol- 
dered into holes in the model skin so that their effective meas- 
uring junctions were located on the inner surface. These 
thermocouples were of No. 36 gage chromel alumel wire and 
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Large heat-transfer model 



'Small heat-transfer model 



Temperature-recovery model 


(a) Heat-transfer and pressure models. 

(b) Temperature-recovery model with 2-inch diameter 
and 0.005-inch-thick wall. 

Figure 2. — Hemisphere-cylinder models. 

were arranged in a spiral of one revolution on the hemisphere 
and of two revolutions on the cylinder. Ten of these thermo- 
couples were located on the hemisphere, fourteen were on the 
cylinder, and one was on the hemisphere-cylinder juncture, 
as shown in table I. 

The pressure-distribution model was 3.000 inches in diam- 
eter, 10.5 inches in length, and 0.10 inch in wall thickness 
and contained 16 pressure orifices located as shown in table I. 
Two of the orifices were located on the base of the model. 
The orifices were 0.040 inch in diameter and were arranged 
in a spiral of one revolution about the sphere in order to 
reduce any orifice-interference effect. The orifices on the 
cylinder were also arranged in a spiral of one revolution. 

The temperature-recovery model constructed of electro- 
plated nickel was 2.000 inches in diameter, 5.00 inches in 
length, and 0.005 inch thick. Temperature measurements 
were taken from 14 chromel alumel thermocouples of No. 30 
gage wire. 

TEST CONDITIONS 

During the tests the tunnel was operated through a 
stagnation-temperature range from 1,040° R to 1,160° R 


TABLE I— ORIFICE AND THERMOCOUPLE LOCATIONS 


Model 

Diameter, 

in. 

Hemisphere 

Cylinder 

Station 

0, deg 

Station 

X, in. 



1 

0 

6 

0. 585 



2 

30 

7 

1. 755 

Heat-transfer 

1. 170 

3 

60 

8 

2. 345 



4 

75 





5 

90 





1 

0 

12 

0. 131 



2 

5 

13 

. 50 



3 

10 

14 

. 75 



4 

20 

15 

1. 00 



5 

30 

16 

1. 50 



6 

45 

17 

2. 00 

Heat-transfer 

3. 025 

7 

60 

18 

3. 00 


8 

70 

19 

4. 00 



9 

80 

20 

5. 00 



10 

85 

21 

6. 00 



11 

90 

22 

7. 00 





23 

7. 50 





24 

8. 00 





25 

8. 50 



1 

0 

8 

0. 393 



2 

15 

9 

. 785 * 

Pressure-dis- 


3 

30 

10 

1. 50 

tribution 

3. 000 

4 

45 

11 

3. 00 


5 

60 

12 

6. 00 



6 

75 

13 

7. 50 



7 

90 

14 

8. 50 

Temperature- 

recovery 

2. 000 

1 

2 

3 

4 

5 

6 

7 

8 
9 

0 

10 

20 

30 

45 

60 

70 

80 

90 

10 

11 

12 

13 

14 

0. 50 
1. 00 

1. 50 

2. 00 
2. 50 


and through a stagnation-pressure range from 10 to 31 atmos- 
pheres. The free-stream Reynolds numbers (based on body 
diameter) were in the range from approximately 0.14X10 6 
to 1.06X10 6 . The calibrated tunnel Mach number at the 
model nose was used to correlate the data. In these tests 
the model nose was approximately 6 inches ahead of the 
center of the test section. At this point the tunnel calibra- 
tion gives a Mach number of 6.82 for a stagnation pressure 
of 31 atmospheres, which agrees with values obtained from 
pressures at the nose of the pressure model. As the pressure 
decreases, the Mach number also decreases slightly, so that, 
at 10 atmospheres stagnation pressure, the Mach number 6 
inches ahead of the test-section center line is 6.74. 

INSTBU MENTATION 

The temperature measurements included stagnation tem- 
perature and model skin temperature. The stagnation tem- 
perature was measured by shielded thermocouples distributed 
at various stations in the settling chamber. 

The potentiometers used for recording the model tem- 
peratures having the greatest transient rate of change were 
cyclic printers with a 1 -second interval and with a 1 -second 
full-scale response. In order to reduce uncertainties as- 
sociated with the time response of the recording instru- 
ments, the thermocouples on the hemisphere were connected 
to be read in order of decreasing temperature about the 
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hemisphere from the stagnation point to the 90° station of 
the hemisphere. The temperature at the stagnation point 
was recorded twice to allow time for balancing out the 
large temperature change between the 90° station and this 
station. Four instruments were connected to read the tem- 
peratures of the thermocouples along the hemisphere in this 
manner and were synchronized so that each thermocouple 
was read once every 3 seconds. 

Model surface pressures were recorded on film by the 
evacuated capsule instruments described in reference- 8. 
The motion of a diaphragm rotates a small mirror to displace 
the trace of a light beam falling on a moving film. Pressure 
cells were chosen to give as near full-scale deflection as 
possible for the measuring station. The stagnation pressure 
was measured with high-accuracy Bourdon tube gages. 

SCHLIEREN SYSTEM 

Two schlieren systems were used in this investigation. 
One system had a single-pass, vertical Z-light path with a 
horizontal knife edge, and the other, a double-pass hori- 
zontal light path with the light source and cut-off knife 
edge at the center of curvature of the concave spherical 
mirror. The double-pass system was used only when greater 
sensitivity was needed for boundary-layer study. Film 
exposures were of approximately 3 microseconds duration. 


ANALYSIS 


HEAT TRANSFER TO MODEL 


The local rate of heat flow into the model was calculated 
from the local temperatures as a function of time and position. 
The general heat-conduction equation (ref. 10) in a homo- 
geneous material is 


dT_k^ 

dr p m C m 


V 2 T 


( 1 ) 


In spherical coordinates this equation becomes: 


dT_ 

r 

11/ 


-F 1 b t 

f ■ Q dT\ 

1 d 2 T"l 

dr 

PmCm L 

r 2 dr' 

V dr) 

1 r 2 sin 6 dd' 

sin u N A j- 
V de / 

r 2 sin d d£ 2 J 


( 2 ) 

Since the model shell is relatively thin witli respect to its 
other dimensions, it is assumed that the first and second 
derivatives of T with respect to 0 and £, as well as the first 
derivative of T with respect to r, are independent of r. 
Equation (2) may then be integrated once with respect to 

bT 

r, the boundary condition — =0 being used at the inner 

surface where r=r in . Then this relation may be solved for 
dT 

-r- to obtain the following: 
dr & 

dT_/ r*-r in * \ Pm c m dT_ 

dr \ 3 r 2 / K dr 



Equation (3) may then be evaluated at the outer surface 
where r=r s and at zero angle of attack so that T is inde- 
pendent of £; this evaluation gives the following: 


fr s *-r in *\ p m c m dT ( r 9 —r in \ d / . 5T\ 

\ dr ) 3 \ 3 r 2 ) k m dr \r s 2 sin d) dd \ 111 ^ dd ) 


(4) 


However at the surface, 



Therefore, 


Qs Pnfim 




—r 




sin 0/ dd 


(si 


sin 0 


dT 

dd 



An examination of the above equation discloses that the 
first term on the right represents q ( , the local net heat 
storage which would be realized in a shell element if the con- 
duction were zero. The term (— ~ 2 ) is approximately 

(Ta—r in ), or the shell thickness, so that in the limit as t—> 0, 
the q t term becomes the same as that used for a flat plate. 
The second term on the right represents the local net heat 
conducted along the model shell. 

When this equation is used to calculate q s at the various 
• • dT 

times and stations, the values of — for the heat-storage 


term may be taken from faired curves of temperature plotted 


against time. 


The values of — 
dd 



for the conduc- 


tion term can be found graphically by taking slopes of faired 
curves. 

In order to give a comparison of the relative size of the 
heat-storage and conduction terms for the tests in which the 


3.025-inch-diameter heat-transfer model was used, typical 
q s and q t values are plotted in figure 3 for a time early in the 
test and for a time late in the test. The difference between 
these curves q s — q t is the value of the conduction term. It 
is seen that, at 5 seconds from the start of the test, the con- 
duction term is only 20 percent of q s for the station at the 
stagnation point; but, near the end of the test, the conduction 
term is about 70 percent of q s at the same station. 

The local heat-transfer coefficients may be determined 
from the standard heat-flow equation for convection, 
q s =h(T s —T r ). The quantities T r and h were assumed 
constant for a given set of test conditions, and h was deter- 
mined from the slope of the straight-line variation of q s 
plotted against T s as is shown in figure 4. In the present 
tests there were some variations in the test conditions during 
a run. The stagnation temperature sometimes varied 20°, 
and the pressure varied nearly 5 percent in the lowest pres- 
sure test. In order to correct partially for the variation in 
stagnation temperature, the values of qJT t were faired 
through a calculated value of T r IT t . The recovery temper- 
ature ratios used in the fairing for the determination of h were 
obtained from the recovery factors for a flat plate given by 
Van Driest (ref. 11) by using the local flow conditions at 
the outside edge of the boundary layer. 


ERRORS 

The present heat-transfer tests were affected by four small 
sources of error which were of uncertain magnitude and were 
not directly incorporated in the computations. One source 
was the heat loss from the model by radiation, and another 
was the heat loss by conduction down the thermocouple 
wire. The third was the assumption that the difference in 
temperature distribution between the outside and the inside 
surfaces of the model was negligible. The fourth source of 
error is the assumption that the heat-transfer coefficient is 
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Figure 3— Heat transfer to the hemisphere of the 3.025-inch-diameter heat-transfer model both with and without the conduction correction. 


independent of the temperature level. The error associated 
with radiation was estimated to be about 3 percent of the 
value of q at the highest temperature condition. The error 
associated with heat flow down the thermocouple wire was 
not estimated directly but was assumed to be small, especially 
on the large heat-transfer model where the thickness of the 
model was 0.10 inch and the diameter of the thermocouple 
wire was 0.005 inch. The temperature difference between the 
outside and the inside skin suface was less than 10° F at the 
stagnation point immediately after the start of a run at j?,=31 
atmospheres. This temperature difference represented an 
error of less than 0.01 in the determination of TJT t . It 


was noticed, however, that each of the first three errors tended 
to reduce the calculated value of the heat-transfer coefficient, 
so that the final results of this data analysis might tend to be 
slightly low when compared with theory. For the assump- 
tion that the heat-transfer coefficient is independent of the 
temperature level, the maximum variation in wall tempera- 
ture occurred at the stagnation point and here the wall 
temperature varied from 180° F to 550° F. By the theory 
of Reshotko and Cohen (ref. 12) a variation in wall tempera- 
ture of this magnitude would result in a deviation of only 
±1 % percent from the mean value of N Nu /^R at the stagna- 
tion point. 
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Figure 4.— Heat-transfer rate plotted against surface temperature for various angles of 6 on the 3.025-inch-diameter heat-transfer model. 
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RESULTS AND DISCUSSION 

PRESSURE DISTRIBUTION 

The longitudinal pressure distributions for the hemisphere- 
cylinder at zero lift are presented in figure 5 for free-stream 
Reynolds numbers of 0.52X10 6 and 1.03X10 6 based on the 
body diameter. These pressure distributions are presented 
as a ratio of local pressure coefficient to the pressure coeffi- 
cient at the stagnation point. This pressure distribution has 
been compared with the modified Newtonian theory from 
reference 13, the data of Stine and Wanlass (ref. 3) at 
M= 1.97 and 3.80, and the data of Oliver (ref. 14) at 
M= 5.8. It will be noted that, at Mach numbers above 3.80, the 
pressure-coefficient distribution curves asymptotically ap- 
proach a single curve, and the local Mach number distribution 
on the hemisphere becomes independent of the free-stream 
Mach number. This was pointed out in reference 3 and is in 
accordance with the hypersonic similarity parameter (ref. 15), 
which predicts that the flow pattern about an arbitrary body 
tends to become stationary as the Mach number increases; 
this tends to occur at a lower Mach number for a blunt body 
than for a slender body. 

In this report the isentropic-flow relations between the 
pressures measured at the stagnation point and the local sur- 
face pressures measured at subsequent stations were used to 
calculate local flow conditions just outside the boundary 
layer. The validity of this method was investigated by the 
use of pitot surveys of the boundary layer. In figure 6 the 
pitot profiles are presented for the hemisphere-cylinder junc- 
ture, a station midway along the cylinder, and a station 
% inch in front of the base of the cylinder. These profiles are 
characteristic of laminar profiles and indicate that the bound- 
ary layer remains laminar over the whole body. The point 
where the pitot profiles break away from the characteristic 
boundary-layer profile corresponds closely to the pitot ratio 
obtained by using the static surface pressure from the pres- 
sure model and the isentropic-flow relations (ref. 16) from the 
stagnation point to the survey stations. Thus, the pressures 
from the body surfaces are sufficient through the use of 
isentropic-flow relations to obtain local flow conditions out- 

22 

side the boundary layer for stations at least up to ^=4.8. 

The local Mach number and temperature ratio as computed 
in this manner are presented in figures 7 and 8. It will be 
noted that the curves exhibit a slight change for the two 
different Reynolds numbers. This was due to a change in 
calibrated tunnel Mach number at different pressure levels. 

LAMINAR RECOVERY FACTOR 

The calculated values of the ratio T T /T t and of the laminar 
recovery factor are shown in figure 9. The recovery factors 
were calculated by a method of Van Driest (ref. 11) and by 
n^jNpr' i- These values are compared with data from the 
temperature-recovery model. Included are the ratios of 
T r j T i and the recovery factors calculated from the uncor- 
rected recovery temperatures or equilibrium temperatures, 
and also from the recovery temperatures corrected for con- 
duction and radiation. 


The laminar recovery factor was also calculated from 
Vr =^J]sjp r) where the Prandtl number was calculated for the 
average of the total and the local stream temperature as was 
suggested by Stine and Wanlass (ref. 3). The results of this 
calculation agreed closely with those from the method of 
Van Driest and were not included in figure 9. 

The experimental values of T r IT t were corrected for con- 
duction and radiation by first calculating the local heat flow 
into the model at equilibrium temperature and then by using 
this value of q s to obtain a value of q s /T t to plot on the 
respective curves of q s /T t against T/T t . In this way a 
point was obtained sufficiently near the T/T t axis to enable the 
fairing from the heat-transfer-model data to q 8 /T t = 0. The 
value of TrlT t obtained by this method was the recovery 
temperature corrected for radiation and conduction. 

The local heat flow q s into the model at equilibrium tem- 
peratures was calculated by using the following equation: 

(7) 

The first term on the right represents the heat flow into the 
model as found by equation (6) for a model in equilibrium, 
hT 

that is, ^—=0. The second term on the right represents the 
Or 

heat flow into the model which is lost by radiation. In a 
slightly different notation this relation is shown in reference 
17. The temperature of the surroundings T a was assumed 
to be constant at 140° F over the entire hemisphere. 

In order to apply equation (7), an estimate of the emis- 
sivity must be made. This estimate depends upon the sur- 
face condition of the material as well as the kind of material 
and its temperature. At the start of the testing the model 
had a polished surface, but after the testing the surface of the 
hemisphere had become slightly roughened. This roughness 
was greatest at the stagnation point and was slight at the 
sphere-cylinder juncture. Because of the uncertainty in the 
value of the emissivity at the elevated temperatures of these 
tests, it was assumed to be constant (0.15) over all the model. 
It was felt that this estimated value of emissivity was not in 
error by more than ±0.05. 

When the corrected values of temperature ratio T r \Tt are 
compared with the values obtained from the two theories, the 
experimental data does not decisively select either theory. 
Since the Van Driest method considers a local Mach number 
effect on the recovery temperature ratio, the values of T r IT t 
obtained by the method were used in the reduction of the 
heat-transfer data. 

TEMPERATURE DISTRIBUTION 

Surf ace- temperature distributions obtained for various 
times during a typical run on the large heat-transfer model 
are presented in figure 10. At the nose the average value of 
the temperature was about 700° R and on the cylinder about 
600° R giving values of TJT m of approximately 7 and 6, 
respectively. 


Pressure coefficient ratio, 
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b 

Figure 5. — Longitudinal pressure-coefficient distribution. 
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Figure 6. — Pitot pressure surveys of flow near the surface of the cylinder. 


Local Mach number, 
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Figure 7. — Variation of local Mach number (outside the boundary layer) along the surface. 
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Figure 8. — Variation of the local temperature ratio (outside the boundary layer) along the surface. 


Recovery factor, rj r , and temperature ratio, 
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Figure 9. — Variation of the surface temperature ratio and recovery factor along the surface. 


INVESTIGATION OF LAMINAR HEAT TRANSFER TO A HEMISPHERE CYLINDER AT M=6.8 


13 



Figure 10. — Typical surface temperature distribution on 3.025-inch-diameter heat-transfer model at various times during a test. 
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HEAT-TRANSFER PARAMETERS N n J^R 

The local laminar heat-transfer parameter was calculated 
on the basis of two different characteristic lengths. For 
one parameter N NlltX l^pR Xi the characteristic length was the 
distance x along the body generator from the stagnation 
point. The air-flow properties used in this parameter were 
the local flow conditions just outside the boundary layer. 
The values of this parameter for M= 6.8 are shown in figure 
1 1 along with the range of experimental data at M= 1 .97 from 
reference 3. These experimental values are compared with 
theories by Stine and Wanlass (ref. 3), by Chapman and 
Rubesin (ref. 18), with a modification of a theory by Sibulkin 
(ref. 7), and with a theory by Reshotko and Cohen (ref. 12). 
For the other heat-transfer parameter N NUtD>a /^R Dt<T the 
diameter of the model was used as the characteristic length, 
and the air-flow properties used were taken as those im- 
mediately behind the normal shock in the center of the bow 


wave. The values of this parameter for M= 6.8 as well as 
the range of values for M= 1 .97 (ref. 3) are shown in figure 12. 
These values are compared with the theory by Stine and 
Wanlass and with a second modification of the theory by 
Sibulkin. 

The local heat-transfer parameter N NUfX /^R~ X) as calculated 
by the theory of Stine and Wanlass, is based upon the local 
Mach number distribution at the outer edge of the boundary 
layer and upon the local Prandtl number. In the present 
tests the Prandtl number varied from about 0.68 at the 
stagnation point to about 0.74 near the end of the cylinder. 
A value of iVp r =0.70 was used throughout the calculation 
of the theory of Stine and Wanlass. This caused an esti- 
mated maximum deviation from the mean heat-transfer 
parameter of only 2 percent. 

In the Stine and Wanlass theory the assumption is made 
that the difference between the actual wall temperature and 
the insulated wall temperature is small. In the present tests 
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Figure 11. Variation of the local heat-transfer parameter along the surface. 
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Figure 12. — Variation of the heat-transfer parameter (based on diameter and conditions behind the normal shock) along the surface. 


this temperature difference varied from about 100° F to 
about 500° F. An examination of figure 11 discloses that the 
data points have a tendency to fall below the theory. This 
difference may be due in part to the rather large tempera- 
ture difference. The reasonable agreement of the data with 
the theory, however, shows that the theory of Stine and 
Wanlass may be used with good results even when the insu- 
lated wall temperature and the actual wall temperature are 
quite different. 

In comparing the theoretical results of Stine and Wanlass 
at M= 6.8 and M= 1.97, a slight discrepancy appears in 
that a discontinuity occurred at the hemisphere-cylinder 
juncture (x/6 = 1 .00) for M= 6.8 and not at M= 1.97. this 
discontinuity can be traced to a difference in local Mach 
number distributions obtained from the two investigations. 
In this investigation the pressure distribution has an abrupt 
change in slope at the juncture (from a high pressure gradient 
to a very low pressure gradient) which gives a similar change 
in the local Mach number gradient (fig. 7). In reference 3 
no such change was noted in the pressure distribution, and 
thus the local Mach number gradient had a gradual change 
at the juncture resulting in a continuous curve for the local 


heat-transfer parameter based on x. Theoretically, a dis- 
continuous pressure gradient is predicted at the juncture. 

Although the use of the theory of Stine and Wanlass leads 
to the calculation of the local heat- transfer parameter 
N NU ' X /iJRz, the parameter N NUtD ,*/^ R d,* may be calculated 
from the parameter based on x, when the free-stream Mach 
number, the local Mach number on the model, and the 
stagnation temperature are known. When the experimental 
and theoretical values of the para meter N Nu>x /^R x in refer- 
ence 3 were converted to N NUtD ^/^R Dtffy a stagnation tem- 
perature of 520° R was assumed. 

The theory of Chapman and Rubesin (ref. 18) is a method 
of obtaining the values of the local heat-transfer parameter 
for a flat plate. The air-flow properties used were those at 
the outer edge of the boundary layer. The distance along 
the body generator was taken as the flat-plate distance from 
the leading edge. This theory depends directly upon the 
air temperature at the outer edge of the boundary layer and 
upon the surface temperature and is not directly a function 
of the Mach number. The theory evaluation at M = 1.97 
is that shown in reference 3. 
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In reference 7 Sibulkin suggested applying his incom- 
pressible solution for stations near the stagnation point of 
a sphere to the nose of a blunt body moving at supersonic 
speeds by using as free-stream conditions those behind the 
center of the bow wave. Since the velocity gradient at the 
nose would not be expected to be the same as in incom- 
pressible flow, Sibulkin’s equation should be used before the 
incompressible velocity gradient around the sphere is 
incorporated, which gives 

Nxu,L = 0.763N Pr OA L 0J' 5 (8) 

The relation was derived with the assumption of incom- 
pressible flow, but the author suggests that the theory may 
be valid for the flow near the stagnation point of a blunt 
body in a supersonic stream. 

Equation (8) may be written as 

hL 

N Nu .l k 

=0.763iV/> r 0 - 4 (9) 

When the characteristic length is x and the flow velocity 
and other physical conditions of the air are evaluated at the 
outer edge of the boundary layer, the first modification of 
Sibulkin’s theory is obtained. This appears as follows: 

hx 

N Nu>x ki 

yiT x / p mx 

V Mz 

=0.763 N Prtl 0A (10) 

The theoretical heat-transfer parameter evaluated in this 
manner is a function only of the Prandtl number. 

When the characteristic length is D and the flow velocity 
and other physical properties of the air are evaluated im- 
mediately behind the normal shock, equation (9) becomes 

hD 

^Nu, D,a ^ 

A H D,a I 

-0.70 (11) 

This equation is the second modification of Sibulkin’s 
theory used in this report. The theoretical heat-transfer 
parameter evaluated in this manner is a function of the 

Prandtl number and of — • 

u a 

In figure 13 (a) the values of ujua are plotted against 
x/D on the hemisphere. These curves show that P is con- 
stant near the stagnation point at a given Mach number. 
In figure 13 (b) experimental values of pD/u * from several 
references, as well as curves computed by using modified 
Newtonian and incompressible flow theories, are plotted 


against Mach number. This curve shows that pD/u ff is a 
function of Mach number. Therefore, N NUtDt ^1l^ 99 as 
found by Sibulkin’s theory near the stagnation point, is a 
function of Mach number. 

When the Newtonian theory was used to find pD/u „ the 
Newtonian pressure coefficient as shown in figure 5 was 
used to determine the pressure distribution around the nose 
of the hemisphere. The correct pressure coefficient at the 
stagnation point to be used with each Mach number was 
found from normal -shock-flow relations (ref. 16). The 
isentropic-flow relations between the pressure at the stagna- 
tion point and the surface pressures computed for subsequent 
stations were used to calculate the local flow conditions just 
outside the boundary layer. This is the same procedure 
that was used to calculate the local flow conditions from the 
measured pressures. 

In figure 11 it is seen that the data indicates a value at 
the stagnation point about 10 percent below the value of 
N Nu ,xl>lR x as found by Sibulkin’s method. In this figure 
it is seen that the data from tests at two different Mach 
numbers, M=1.97 and M= 6.8, agree within the scatter of 
the data. In figure 12 it is seen that the data at the free- 
stream Mach numbers of 1.97 and 6.8 indicate values which, 
within t he sc atter, are consistent with the predictions of 
N NUt D ,J^R D ,a by Sibulkin’s method. It is important to 
note, however, that the variation of pD/u 9 with the free- 
stream Mach number must be taken into account and that 
the prediction of P by the potential theory may not be used 
for a hemispherical nose subjected to supersonic flow. 

LOCAL STANTON NUMBER 

In figure 14 the local Stanton number is plotted against 
the local Reynolds number. The data presented includes 
points from 0= 5° to 0=45° on the hemisphere surface. 

By assuming that the velocity and temperature boundary- 
layer profiles were similar and of parabolic shape, Sibulkin 
(presented in ref. 2) derived an approximate solution for 
N NUiD over the entire hemisphere for incompressible flow 
which, when transformed for the local conditions outside 
the boundary layer and distance from the stagnation point, 
becomes 

N„ u . x =k(d)x(y'*N Pr . ,«•■> ( 12 ) 

where k(6) is a proportionality function which depends on 
the angular position of the point under consideration 
measured from the axis of flow symmetry. 

If P is assumed constant (a valid assumption for 0 up to 
45° in these tests, see fig. 13), the equation reduces to 

N Nu , x =k($)R x 0 - 5 N Prt t 0A (13) 

From reference 6 it may be seen that k{6) may be approxi- 
mated near the stagnation point by the value of 0.763, and 
the local Stanton number near the stagnation point becomes 


N: 


St, l 


N Nu,x 

RxN PTt i 

=0.763/i > z -0 * W/v t i 


- 0.6 


(14) 


INVESTIGATION OF LAMINAR HEAT TRANSFER TO A HEMISPHERE CYLINDER AT M = 0.8 


17 

























































































. 


— 

— 



— 

— 


, 




— 



— 

— 











r 





















— 




> 











==- 































































< 

> 

O 









































A 





































□ / 















































































A 





































o 

/ 








































/A 







































/ 






























"'ln< 

:om 

pressible theory 

7 

a 































/ 








r~ 

































' \ 

> 

\ 





/ 


































a 


\ 



/ 

*■ — 
















o Stine and Wanlass 

□ Stalder and Nielsen 

O Korobkm 

a Korobkm 

v Present tests 

1 7 Oliver 









s~ 

a 

\ 

\ 

L 

v 

























7 























r 
















































M 

odified Newtonian theory 

IVI 
























Kayieign- Janzen meinoa, rer.io 



















































































































































































_L 






















(b) 

c 


2 3 4 5 6 

Free-stream Mach number 


(a) Nondimensional local flow velocity along the surface. 

(b) Nondimensional velocity gradient at the stagnation point of the hemisphere plotted against free-stream Mach number. 
Figure 13. — Nondimensional flow parameter on the hemisphere outside the boundary layer. 
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Local Reynolds number 

Figure 14. — Variation of the local Stanton number with the local Reynolds number for values of 0 from 0° to 45°. 


For this equation to give the correct Stanton number around 
the hemisphere surface, dujdx or 0 must be constant and 
k(d) should not depart radically from the value of 0.763. 

In figure 14 the heat-transfer data has been compared 
with Sibulkin’s stagnation-point solution shown in equation 
(14) . As seen from figure 13 (a) , can be considered constant 


up to values of 0 greater than 45° (j^=0.4^- In figure 1 1 


an 


indication of the variation of k{6) may be obtained from the 
decrease of the heat-transfer parameter from the constant 
stagnation-point value as x is increased. Indeed, even the 
difference in level which is observed between the stagnation- 
point curve and the data in figure 14 may also be detected 
in figure 11. Sibulkin’s stagnation-point value is approxi- 
mately 9 percent above the data at the nose and about 15 
percent above the data at the 45° station. Above 45°, k(d) 
diverges rapidly and 0 begins to diverge. Inasmuch as the 
data above 45° do not correlate, they are not presented in 
figure 14. 


AVERAGE HEAT-TRANSFER COEFFICIENT 

The average heat-transfer coefficient determined experi- 
mentally in this investigation and by Staid er and Nielsen in 
reference 4 are presented in figure 15. The data from Stalder 
and Nielsen are those for M=0.12 to 0.17, M= 1.75, 
M= 2.67, and M= 5.04. Also included in figure 15 are exper- 
imental points found by an integration of experimental data 
by Stine and Wanlass (ref. 3) at M= 1.97. The data has 
been compared with values found by using the theory of 
Stine and Wanlass at M=1.97 and M= 6.8 and with the in- 
compressible theory of Sibulkin used by Korobkin (ref. 2). 
Korobkin adjusted the latter solution to pass through the 
stagnation point for incompressible flow (ref. 7). 

The data from Stalder and Nielsen (ref. 4) for Mach 
numbers below 2.67 appear to correlate with the adjusted 
incompressible theory. The data of reference 4 were evalu- 
ated, however, in a different manner from the data of this 
report. In the present report the local values of the heat- 
transfer coefficient were calculated from the local heat con- 
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Reynolds number behind normal shock 

Figure 15. — Average heat transfer on the hemisphere plotted against Reynolds number for conditions behind a normal shock. 


vection and the local recovery temperature. Then these 
local values of h were integrated over the area of the hemi- 
sphere to find the average value. In reference 4 the average 
heat-transfer coefficients were evaluated by measuring the 
total heat convection to the hemisphere and an “effective” 
average recovery temperature found by an extrapolation to 
zero net heat transfer to the surface. The different results 
as found by the two methods are seen in figure 15. 

The integrated data from Stine and Wanlass and the 
data from this investigation show a good correlation at 
Mach numbers of 1.97 and 6.80 and Reynolds numbers 
from 0.02X10 6 to 1.5X10 6 . 

SCHLIEREN OBSERVATIONS 

The normal shock and body boundary layer are shown in 
the schlieren photographs (fig. 16) for a single-pass and 
double-pass schlieren system. The relatively nonsensitive 
single-pass system clearly shows the normal shock while 
the sensitive double-pass system was used to show that the 
boundary layer is laminar on the surface of the cylinder. 
The two flow photographs appear to show additional dis- 
turbances about three-fifths of the distance along the cylinder 
from the hemisphere-cylinder juncture. These disturbances 
result from the intersection of the bow shock wave and the 
boundar}^ layer at the windows. The heating effect was 
great enough so that a similar picture of this region was 
obtained immediately after a run, as shown in the post-run 
zero picture for the double-pass system. 


CONCLUSIONS 

A program to investigate the aerodynamic heat transfer 
to a nonisothermal hemisphere-cylinder has been conducted 
in the Langley 11-inch hypersonic tunnel at a Mach number 
of 6.8 and a Reynolds number based on diameter and free- 
stream conditions from approximately 0.14X 10 6 to 1.06X 10 6 . 
In this investigation the boundary layer was laminar having 
average values of T s /T m at the nose of about 7 and about 6 
on the cylinder (where T s is the local free-stream temperature 
just outside the boundary layer and T m is the free-stream 
temperature ahead of the normal shock). The data have 
been correlated with theoretical analyses and results of other 
experiments at lower Mach numbers. The following results 
were obtained : 

1. The experimental heat-transfer coefficients from this 
investigation on a nonisothermal (constant wall thickness) 
hemisphere-cylinder with heat flowing to the surface were 
slightly less over the whole body than those predicted by 
the theory of Stine and Wanlass (NACA Technical Note 
3344) for an isothermal surface. 

2. A modification of Sibulkin’s stagnation point solution 
gave the trend of the local Stanton number with local 
Reynolds number for angles up to 45° from the stagnation 
point. The calculated values, however, were approximately 
12 percent higher than the experimental values. 

3. Pitot profiles taken at a Mach number of 6.8 on a hemi- 
sphere-cylinder have verified that the local Mach number or 
velocity outside the boundary layer required in the use of 
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Pre-run zero picture Flow picture 

Single-pass schlieren system 



Post-run zero picture Flow picture 

Double-pass schlieren system L-92450 


Figure 16. — Schlieren photographs of shock and boundary layer on hemisphere-cylinder at M = 6.8. 


the theories may be computed from the surface pressures 
by using isentropic flow relations and conditions behind a 
normal shock. The experimental pressure distribution at 
a Mach number of 6.8 is closely predicted by the modified 
Newtonian theory. The velocity gradients calculated by 
using modified Newtonian theory at the stagnation point 
vary with Mach number and are in good agreement with 
those obtained from measured pressures for Mach numbers 
from 1.2 to 6.8. 

4. At the stagnation point a second modification of the 
theory by Sibulkin using the diameter and conditions behind 
the normal shock was in good agreement with the experiment 
when the velocity gradient at the stagnation point appro- 
priate to the free-stream Macli number was used. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 21 , 1956. 
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